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SUMMARY 

A source-distribution method is applied to derive a solution for 
the time-dependent surface velocity potential of thin finite vlngB at 
supersonic speeds. The solution is illustrated by evaluating the 
upwash over the tip of an arbitrary-plan-boundary wing having a 
supersonic and subsonic leading edge. The upvash then gives the 
effective sources of the flow field lying between the wing plan 
boundary and the foremost Mach waves, which are applied in the 
derivation of the wing-surface velocity potential. A simple example, 
the load distribution on a wing whose effective angle of attack is 
changing linearly with time, is Included to illustrate the appli- 
cations of the derived expressions. 


INTRODUCTION 

The analysis of the aerodynamic effects in the vicinity of thin 
wings at supersonic speeds can be simplified to obtain useful results 
by means of the linearized theory. Steady-state or time-independent 
solutions have been obtained for enough cases (for example, refer- 
ences 1 to 27) that the essential features of the load distributions 
on various parts of the wing can be determined by graphical or 
analytical methods. The time-dependent load distributions are more 
difficult to obtain. Such problems include the transient effects 
of gusts, changes in angle of attack, skin vibration, and flutter. 

A number of investigators have studied two-dimensional time- 

dependent flows over thin wings. ' These flows are generally included 
as special cases of the theory of reference 27. The method of 
reference 27 is similar to the steady-state solution of reference 1 
and includes three-dimensional or finite wing solutions in cases 
where the aerodynamic effects of the bottom and top wing surfaces are 
independent. No solutions are known to have been published for 
cases involving interaction between the flow over the bottom and 
top wing surfaces. 


/ 


2 


MCA TN No. 1699 


The concepts of references 9 and 18 that lead to solutions, 
which Include the effects of Interaction between the bottom and 
top wing surfaces in the steady state, can also be applied to 
obtain time-dependent solutions. The derivation of the surface 
velocity potential in regions Influenced by subsonic leading 
edges was completed during January 1948 at the MCA Cleveland 
laboratory and is presented herein. 


AMLTSIS 

In order to unify the discussion, parts of the fundamental 
treatment presented in reference 27 are repeated. The analysis 
includes the derivations of (l) the time-dependent, linearised 
partial-differential equation for the perturbation-velocity 
potential of an ideal fluid, (2) the fundamental solution that 
will satisfy the boundary conditions on the wing, (3) the upvash 
between the wing boundary and the foremost Mach line, (4) the 
velocity potential on the surface of the wing. 

Differential equation. - The linearised Euler's equations 
for a compressible fluid may be written 


d5p 

+ TJ—? m - 

..LSe 

dxdt 

ax 2 

P 0 

d 2 * 

dydt 


-^ly 

d 2 op 

_a 8 i> 

_ _L 

dsdt 


’ Oo Si 


where 

cp perturbation-velocity potential 

t time 

TJ free-stream velooity 

p Q free-stream density 


( 1 ) 



BACA TN Bo. 1699 


3 


p static pressure 

x,y,z Cartesian coordinates (free stream parallel to x axis) 

(A complete list of symbols is Included In appendix A. ) 

When equations(l) are multiplied by dx, dy, and dz respec- 
tively, added, and Integrated, the result is 




dt dx 


( 2 ) 


where g(t) is an Integration constant at any given time, 
The linearized continuity equation can he written 




(3) 


or, because the velocity of sound c is 


equation (3) becomes 



Substitution of p/p Q from equation (2) in equation (3) yields 

( X .^)^e + ^2 + ^2.i^e + 25^SL + i5a (4) 

3x 2 8y 2 o 2 dt 2 o 2 arSt e 2 at 


where M is the Mach number of the free stream and the zero sub- 
script has been dropped from c. 

Equation (4) is the required linearized partial-differential 
equation for the veloolty potential. If qp is independent of time, 
the Prandtl-Glauert equation immediately results. 
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The function g(t) depends on the condition of the flow ahead 
of the body. If the flow Is uniform and undisturbed, g will be 
constant (from, equation (2)) and equal to Pq/Pq* function g 

will be assumed constant In the rest of the analysis. 

A change In variables will convert equation (4) to a standard 
form of the wave equation. The transformation equations are 

x‘ - x 

y* - Jl - M 2 y 
2 * - - V? 2 

t« O (1 - M 2 ) t (5) 

Application of equation (5) In equation (4), with g set equal to 
a constant, gives 


a 2 q> ^ a 2 <p ^ a 2 ? JL*?SL 

dx* 2 Sy* 2 S*« 2 c 2 dt' 2 


( 6 ) 


Basic solutions of equation (6) corresponding to spherical waves 
are 



andqpB f 
r ' 




( 7 ) 


where r* 
function. 


V*' 2 + y' 2 + s' 2 , P 2 


M 2 - 1, and f 


is an arbitrary 


The basic solution for the supersonic case Is obtained as the 

sum of equations (7). (See reference 27.) If this solution Is 
transformed to a general point In the- x,y,s spaoe, the basic 
solution of equation (4) assumes the form 



where £,t) ; and t are Cartesian coordinates In the x,y, and. z directions, respectively. 
Beoause the operations Indicated In equation (4) are Independent of t,i), and f , f nay 
be a function of these variables. Extended solutions of equation (8) may be obtained by 
integration with respect to any of the variables |,t|, or£. 

For thin wings, the souroes of the disturbances will lie In the plane of the wing, which 
nay be defined by the relation £ = 0. The problem Is then to determine the function f In 
terns of known or computed wing or streamline deflections. 


Required, source strength for thin wings. - In order to shorten the equations, the 

quantity V ( x -|) 2 - p 8 (y-n) 2 - P 8 b 2 nay be replaced by the symbol r. Equation (8) then 
becomes 




(8a) 


The porturbatlon-veloclty component dqjto* nay be obtained by differentiation to give 


05 


da r 



( 9 ) 


The vlng dictates that dq^da is pot zero near a ■ 0. Equation (9), however, shove that the con- 
tribution of the sources to- &p/da Is zero, except- for those souroes near r ■ 0. The quantity r 
is nearly aero either for snail values of (x-|), (y-t)), and s or for (x-S) 2 + a 2 3« 

Tbs second condition gives a lover-order aero in the denoolnator of equation (9) than the first con- 
dition and does not contribute to the value of dp/da at the point (x,y,o) on the vlng surface. 


Bear the point (x,y,0) the quantities (x-i) and r are nearly aero. Sis function f then 
reduces to a function of x.y, and t. The veloolty potential at point (x,y) due to the eubelenenbs 
in the vioinity of the point can then he obtained by integration of the souroes Included in the forward 
Jfaoh cone. Fran figure 1, the field of integration is hounded hy the curves | ■ |j, and 

(x- g) 2 - 0 2 (y-tj) 2 - p 2 a 2 m 0. The velocity potential la then 



- - p. -| X ) 


( 10 ) 


Partial, differentation of equation (10) with respect to a gives 


- w(x,y,t) - - 2*f(x,y,t) (11) 
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where w Is the looal z component of the perturbation velocity. 
Equation (ll) defines the strength of the local source in terms of 
the z component of the perturbation velocity. The fact that equa- 
tion (ll) is independent of the value of is in agreement with 
the statement that only those subelements near the point (x,y,0) 
contribute to the vertical velocity dq/Sz on the wing surface. 
Physically, this evaluation infers that the flow is tangent to the 
wing surface at each local point. 

The time-delay terms of equation (8) can be denoted by T ft and 
t ^ where ' * 



The velocity potential at any point (x,y) is obtained by integrating 
the sources in the z ■> 0 plane over the area S included in the 
forward Mach cone. By use of equation (ll), the velocity potential 
becomes 


Equation (13) was derived by Garrick and Bubincw in reference 27 and 
gives the velocity potential at any point in space in terms of the 
z component of the perturbation velocity in the z = 0 plane. The 
equation reduces to the steady- state solution of Puckett (reference l) 
when w is independent of time. 

A physical interpretation may be given for the time delays T & 

and. T ^ of equations (12) and (13). If a disturbance is generated 
at point (£ ,tj) at time t = 0, the wave front from that disturbance 


[w(l,*1,t-T a )+ w ( gyT^t-Tfr)] d|dn 

V(x- 1) 2 - P 2 (y-t|) 2 - e 2 * 2 
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will travel outward as a spherical wave about a center that moves with 
the free- stream velocity. (The trace of the waves on the z ■ 0 plane 
is illustrated in fig. 2.) The wave front will enter and emerge from 
the point (x,y,») at two later times T a and T-j,. T he equation of 

the spherical wave path that passes through the point (x,y,z) is 


(x-£-UT ) 2 + (y-tj) 2 + s 2 « c 2 *! 2 (14) 

The solution for T is 

t ''l(x-t) 2 - p 2 (y-n) 2 - 

P 2 C P 2 C 

which gives equations (12). (At a given point (x,y,z), the strengths 
of the same wave at the two times T a and are equal despite the 

change in the radius of the wave front.) At a given time t, only the 
wave fronts that are entering and emerging from the point (x,y,z) con- 
tribute to the velocity potential. These two waves originated at the 
point (|,t),0) at times (t-T a ) and (t-Tjj) . 

The remainder of the analysis is primarily concerned with the 
aerodynamics in the plane of the wing so that z may be set equal to 
zero. Equations (12) and (13) may also be conveniently expressed in 
an oblique coordinate system whose axes lie parallel to the Mach lines. 
(See fig. 3.) The transformation equations are 



V 


E W 

n 

= i (▼-») 

”» * # (i - p7) 

v w 

- is (l+?y) 

1 <v ra w> • 

7 



(15) 


Inasmuch as the elemental area in the (u,v) coordinate system is 

•2j| dudv. equations (13) and (12), in the case of z * 0, became 
M 2 



MCA TN Ho. 1699 


9 




M( V y +V u) - 2 Y(\-u)(v v -v) 


(17) 


Equation (16) gives the velocity potential in the z « 0 plane 
in terms of the perturbation-velooity component v normal to the 
plans. If only supersonic leading edges are Included in the forward 
Mach cone from (x,y), v may he evaluated in terms of the effective 
-wing slopes O measured in tj * constant planes hy the relation 

v « TSa (18) 

& 

If a subsonic leading (or trailing) edge is also included in the 
forward Mach cone from (x,y), the slopes of the stream-lines X 
associated with the upwash between the wing, boundary and the foremost 
Mach wave must be evaluated and included in the calculation (equa- 
tion (16)) for the velocity potential. 

Upwash between wing boundary and formoat Mach line. - The slopes 
of the streamlines in the region Sp of figure 4 could conceivably 
be generated by a thin wing or diaphragm, as employed in reference 9. 

The flows above and below the z * 0 plane may then be independently 
treated. The velocity potential on the top surface of the z * 0 
plane is given by equations (16) and (18) as 
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<P„ 





duly 


yCup-uJCvD-y) 


_X(u,y,t-T a ) + Xu,y,t-T b )] dpdv (ig) 

y(t^-u)(v D -v) 

y her e up and y D are the coordinates of the point at which q> la 
evaluated, G ^ represents the slopes of the strea m l ine s on the top 
wing surface, arid A represents the slopes of the stre amline s in the 
field Sp from the point of view of the top wing surface. S imi larly, 
the potential on the bottom surface of the z = 0 plane is 




^ B (n,y,t-r a ) + gju,T,t 7 Tj] dudv 
^ (up-uKvjj-y) 



-|x(u,y,t-T a ) + Mn^t-T^) J dudv 
if (ujj-uJCv^v) 


(20) 


The pressure at a given point of the field Sjj (fig. 4) In the 
plane of the wing can be calculated by substituting either cpij> or 
cp B into equation (2). The two computations of the pressure can ■ 
then be equated: 



MCA TN No. 169J9 


11 


Equation (21) has the solution - 

<P T -cp B + 2H(x-Ut,y) (21a) 

where H is an integration function. Substitution of equations (19) 
and (20) in (21a) yields 


-4 ff . - -j- / r *a l^^ +h 

J J V(^>" u) ( v D“ v ) J J s 2 V(%‘u)(v d -t) 

U US D S w 

(21b) 

(In equation (21b) and in sene of the equations to follow, the 
notation !\ a means Mu,v,t-T & ) and X a ^ means 

[x(u,T,t-r a ) + >(u,v,t-T-jj)J. A similar notation is applied to 
Ofe, a T , and (Oq-Oj,).) If equation (21b) is substituted into (19), 
the result is 


*T 



( VVa.b dlldT 

2 'l/C^j-u) (▼]>-▼) 


+ H 


( 22 ) 


Equation (22) represents the Telocity potential in the plane of the 
wing for the region Sp. Similarly, formulation of would change 

only the sign affixed to H. The function 2H of equation (21a) 
represents the difference in potential across the z « 0 plane, 
corresponding to the strength of vortioity in the wake of the wing 
(reference 18) . For flat-plate wings, Og + Cj ■ 0 and H is 

just the potential on the top surface of the vortex sheet. 

The foremost Mach wave (fig. 4 or 5) originating on the leading 
edge generally represents a line of infinitesimal disturbance along 
which H(x-Ut,y) can be set equal to zero at all times. The func- 
tion H remains zero along y « constant lines for values of x 
not intercepted by the wing or a material body (region Sp ^ of 
fig. 5). The region Sp 2 of figure 5 generally contains a vortex 
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sheet lying In the plane of the wing and H Is not zero. The func- 
tion H(y), established along the wing trailing edge at some time t, 
remains unaltered for later 'times along a curve that sweeps down- 
stream with the free- stream velocity and has the form of the wing 
trailing edge. The rest of the discussion is concerned with only the 
effects of leading edges, that Is, those cases for which H is zero. 

The origin of coordinates for the wing shown on figure 4 is 
placed at the junction of the supersonic and subsonic leading edge. 

The supersonic and subsonic leading edges are respectively defined by 
the equations 


▼j(tt) 

or 

U a U-^(v) 


v 2 (u) 

or 

U a U^V) 

(23) 


Equation (21b) with H set equal to zero, then becomes 


P^D 


UO 


du 


^ D' u 


n v D 


Uv 9 (u) 


X >dv 

AfV* 


P%> 


du 




0^2 (u) 


Jvi(u) 


(Q B - CT T) a ,b dY 


(21c) 


In the steady-state solution, X and (Og-CSp) are Independent 
of Up (references 9 and 11 ), so that the Integrations with respect 
to v may be equated to give 


OTt 


pv 2 (u) 


lvp(u) 


X a,b dY 


'■^(u) 


( < ^-q r )a,b dT 
2 y d- y 


(24) 


The reduction of equation (21c) to give equation (24) for the 
time- Independent case shows that only those wing slopes along 
u *s up contribute to the upwash for points on this Mach line. In 

the time-dependent case, both X and (Oj-Crji) contain functions 
of UD, and the validity of the reduction may be questioned. If 
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in the time-dependent case only those wing slopes along u c uj> 

were to contribute to the upvash for points on this Mach line, 

X and ( 023 -C 3 O. would heccme Independent of Uj) and the reduction 
of equation (21c) to equation (24) would again he justified. An 
argument to substantiate this choice is presented in appendix B. 

If u is replaced by uq In equation (24), the two time 
delays become equal and are linear with respect to v: 


T b 


Pc 


Each Inlniteslmal wing element then produces at time t and place Vj> 
an Increment In upwash corresponding to the steady-state effect of 
wing elements whose slopes are the same as the time-dependent elements 
evaluated at time t-T & . The solution of equation (24) is the sum 

of 1 a series of infinitesimal steady- state solutions, each of which 
satisfies Abel's equation (reference 11) and each of which requires 
Integration only over the wing slopes along u = up. 

In order to Illustrate the argument, the Increment in X at 
point (up, v D ) (fig. 6 ) due to a steady-state wing element of 

length dv at point (ujj,v) may be obtained (appendix B) from 
Abel's equation In the manner of reference 11 as 


dX (Pb-Pt) 4 y 2“ t 

dv = 

2 rt(v D -v)/y v D -v 2 


(25) 


Equation (25) may be applied In either the time-dependent or the 
time-independent cases. For the time- Independent cases, the wing 
slopes are evaluated at uq,v. For the time-dependent solutions, 

the wing slopes are evaluated at (un,v) at time t - 

^ pc 

Integration of equation (25) across the wing gives 
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where ,,, means (0^-0^,). That equation (26) Is a solution of equa- 
tion (24) is shown in appendix B. 

Evaluation of velocity potential on wing surface, - Equations (19) 
and (26) now allow the calculations of the velocity potential on the 
wing surface. With reference to figure 7, the velocity potential at 
the point is 



(27) 

For the wing of figure 7, the integration of X over the 
area Sjj is » 



If X from equation (26) is substituted into equation (28), there 
results 



fLG 





(28a) 

\ 

The Integration of X over the flew field Bp can he replaced by the integration of a 
function over the wing area 8*3.. This replacement Is accomplished by Interchanging 
the order of Integration of v and v D . Equation (28a) then becomes 



(28b) , 
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Equation (2Sb) represents the contribution to the velocity potential of the upwash error 
the wing tip and. con he used, to rewrite ports of equation (27). If the -velocity potential 
from equation (27) is written in full for the wing of figure 7, the following results is 
obtained.: 



(continued on the following page) 
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Although equation (27a) loots rather formidable, there are 
actually only two types of t era. Simple variations are obtained, 
by changing signs from minus to plus, by replacing bottca-surfaoe 
wing slopes by top-surface wing slopes, or by altering the limits 
of integration. Simplifications also occur in specific examples, 
although the indicated integrations to obtain explicit solutions 
are generally difficult to perform. 

If the wing slopes In the vicinity of the wing tip vary con- 
tinuously with either time or position, equation (27a) can be 
applied in its present form without altering the limit s of integra- 
tion. The equations for the wing slope would be Inserted into 
equation (27a) and the indicated Integrations conducted. Thi s type 
of calculation has been illustrated for perlodio oscillations of 
flat-plate wings in references 27 and 28 when the bottom- and top- 
wing slopes are independent. A less complicated example, in which 
the angle of attack of all wing elements of a flat plate varies 
linearly with time. Illustrates this type of calculation. 

If the angle of attack of all wing elements is changed at a 
uniform rate m for all times, the wing slopes may be expressed 
as 


CFg o a + mt m — o j 


(29) 


where a is the angle of attack of the wing at time t » 0. 

These effective wing slopes may correspond to a constant accelera- 
tion of if In the z direction. Substitution of equation (29) 
into equation (27a) gives 


(TO W 






du 


UU^V ) 


Vv“ 


dv 


V. (u) 




« 


P«w 


rv*« 


Urn I 
Mtt 

U 


du 


*2< v w> 




liv^u) 


V Bo Be / 


(30) 
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Integration with respect to ▼ produces 


<P, 


2TJ 


T Mrt 


— 

n«w 

P«w 

(a+nrt) 

V T w-» « 


A/VU • 3P ° 


_ L 

”2 > y J 

*2 


(Y v -7l(u)) 

i/v 5- 


4u 


m 

pc 








(30a) 


If the equation for the leading edge of the wing is v = - k^u, 
equation (30a) may he integrated to obtain 


tp_ *= .22 < 

T ifa 




ll^ 1 ' t )"2 V 


(Uy-UgKVy ♦U^) 


(a +mt) (k^Uy +y v ) - 


(ki + lJaUiUy +v w )' 


tan 


-1 


4 P°*L J V T w +k l *2 

(31) 

The pressure coefficient nay he derived from equation (2) as 


, _ p - p o- 

'P 1 2 


u\ 


It) 


(32) 


Substitution of equation (31) into equation (32) results in the pres- 
sure coefficient on a family of wings in the region influenced hy 
the subsonic leading edge. 
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(53) 

where ug(v) Is evaluated at v * v y . (Although equation (33) was 

derived for the case of constant acceleration nftJ In the z direc- 
tion, the solution may he combined with other equations to evaluate 
the pressure distributions for a variety of wing motions. For 
example, a wing rotating with constant rate m about an axis of 
pitch fixed with respect to the wing and lying in the 2*0 plane 
would have the pressure distribution given by equation (33) superposed 
on the pressure distribution associated with uniform rates of pitch 
described in reference 21. ) 

The steady-state solution included in reference 18 is obtained 
from equation (33) by setting m = 0. . The solution far the infinite 
swept wing is obtained by setting v v * u 2 (v w ) * 0 and computing 
Cp. The values of C p obtained along the line v w = 0 are then 
constant along lines parallel to the leading edge. The load dis- 
tributions of a family of wing plan boundaries may be evaluated by 
choosing the desired equation u * ug(v) for the wing-tip plan 
boundary. 

The general equation (27a) also includes solutions for which , 
a finite number of discontinuities can exist in the wing slopes with 
respect to either position or time. The procedure is the as 
in the case of continuous wing slopes, exoept that the fields 


974 



MCA TN No. 1699 


21 


of Integration will "be subdivided In accordance with the requirement 
of the discontinuity. For example. In the problem of determining 
gust loads, the effective angle of attack of the wing will change 
dlsconblnuously along a line parallel to the y axis that moves 
downstream with the free- stream velocity. If the Initial wing slope 
Is taken as zero and the Increment In wing slope due to the gustf Is 
a, the slopes of the wing at time t - T (as time appears In equa- 
tion (27a)) are either zero or a according to the inequalities 

|> U(t-T), * “°T ■ 0 

| <U(t-T), g b c -C T «= a 

The two-dimensional solution of this problem is presented in ref- 
erence 29. The three-dimensional solution of this and other transient 
problems within the scope of equation (27a) requires further research. 


Flight Propulsion Besearoh Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio , May 18, 1948. 
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APPENDIX A 
SYMBOLS 

She following symbols are used In this report, 
pressure coefficient, 

velocity of sound 
function 

integration function of time (herein considered as 
constant ) 

(integration function of x-Ut and y 
Integral 

constant greater than zero 
free- stream Mach, number 

i 

rate of Increase of wing angle of attack 
static pressure 

V (*-£) 2 - P 2 (y- 1 i ) 2 - P 2 z 2 

. ^x*2 + y* 2 + z' 2 <* 

plan-form area 
manipulation variable 
time 

(l-H^ t + f 

free- stream velocity, taken parallel to x axis 

oblique coordinates whose axes lie parallel to Mack 
lines in z » 0 plane 


P“Pp 
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w z component of perturbation Telocity, (taken as 

positive in direction of an outwardly drawn normal 
from plane of wing) 

x.y,z Cartesian coordinates 


x'jy'jZ' transformed Carte sian coordinates 

(x’ = x, y* ssVl-M^y, z’ «*/l-M^z) 

|,t),£ Cartesian coordinates in x,y,z directions, respectively 


a 

0 


angle of attack 


cotagent of Mach angle. 



1 


\ effective slopes of streamlines (measured In T] * constant 

planes) in z = 0 plane between wing boundary and 

foremost Mach line, A ^ — 

7 TT 


p density 

O effective wing-section slopes measured in tj = constant 

planes ( d = w/U) 

T time delay 

cp perturbation-velocity potential 


Subscripts: 


0 

1,2,3 . 

l 

a,b 

B 

D 

T 


free stream 

numbered areas or wing-plan-boundary equations 
time delays T & and T-jj 

bottom (of wing) 
upwash. field 
top (of wing) 


w 


wing 
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Examples: 



V 1 

curve v = v^(u) 

along supersonic leading edge 

«2 

curve u s ^(v) 

along subsonic leading edge 

S w(l + 2) 

wing areas 1 plus 

2 

&Ug 

dv 

w 

derivative of curve ^(v^) with respect to 

^a,b 

slope at time t- 

T & plus slope at time t - T )> 


( 03 — Oip) a ^ difference "between "bottom and top wing slopes at 
* time t -T a plus this difference at time t- T -j, 

°B,T °B ~ °T 

U^/v oblique coordinates of point x,y on wing 



MCA TN No. 1699 


25 


APPENDIX B- 
EVALUATION OF UPWASH 

Origin of upwash. - The fundamental source solution (equa- 
tion (8); for the three-dimensional velocity potential of thin 
wings moving at supersonic speeds can he written as 

cp = . (Bl) 

T 

If this equation is partially differentiated with respect to z, 
the result is 

*P = 22 S’? (B8) 

dz Sr dz r dr 

If z is allowed to approach zero, the quantity Scp/Sz will 
approach zero unless ± djp approaches infinity (that is, r 

r Sr 

approaches zero) . The fundamental solution (Bl) can then he Integrated 
over a surface of sources in the z » 0 plane to obtain an extended 
velocity potential. Because of the form of equation (B2), the z com- 
ponent of the perturbation velocity generated by this potential will 
arise only from those sub elements near r «* 0. 

On the wing, the flow must be tangent at each point to a defined 
surface. The quantity Sq/Sz is therefore determined at each local 
point by the wing. Under such circumstances and because of the wing 
restraint, only the subelements near r *« 0 obtained by assuming 
(x-!)ss(y-Ti)=s(z-£) contribute to SfiP/Sz . No such restraint 
exists for the upwash field. If the same condition is imposed to 
evaluate dq/dz in the upwash field part of the z *= 0 plane, the 
trivial answer _ S*£ results, 
dz dz 

A fundamental distinction exists between the evaluation of 
dq/dz on the surface of the wing and in the upwash field. The 
wing surface is a restraint and can generate primary impulses that 
lead to the velocity potential and hence to the flow distributions. 

The upwash field is unrestrained, however, and can transfer from 
the bottom to the top wing surfaces only those impulses that have 
already been generated by the wing. 
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The upwash most arise, however, from, those sub elements near 
r ■ 0. Because the unrestrained subelements In the vioinity of the 
point (x-|) * (y-rj) s z m 0 cannot generate upwash, dcp/dz 
(for z s 0 ) must he generated hy those subelements In the vicinity 

of the curve - p^(y-q)^ = 6^z^ * 0. This condition Implies 

that only those wing slopes along u = Up generate the upwash for 
points on this Mach line. The defining equation for X should then 
he restricted hy setting u = Uq, for either time- independent or 
time-dependent flows. This result was proved directly from the 
Integral equation (21c) for time-independent cases (reference 11). 


Differential equation for X. - The Increment dX at point 
Uq,Vjj due to the wing elements dv (fig. 6) of wing slopes 

(Xt> _Qm 

— s — may he determined from equation (6) of reference 11 as 


0® 


dX«! JL 
X os 


dTr 


V s 


nv+dv 


(Cg - a T )dv 


(B3) 


2a^v^ - v 


where s Is the manipulation variable and Is evaluated at s > v_. 
The integration of dv over an infinitesimal distance results ° 
In the removal of one Integration sign. Also, (q B - a T ) may he 
considered as constant. Equation (B3) then becomes 


dX 

dv 


a B' 


2x 


^8 






dv T 


^(s - v D )(v D - v) 


(qg - qip) - ▼ 

2 it (s - v)^a - v 2 


(B3a) 


Replacement of the manipulation variable s hy v B gives equa- 
tion (25) 



( g B ~ 
2* 


V ? 2^ r 

(v D - v)Yv d - V 2 


(B5b) 



974 


MCA TN No. 1699 


27 


Demonstration by substitution that equation (26) satisfies 
equation (24) when u ■ Up. - if the Integration variable v Is 

changed to s to prevent confusion, equation (26) assumes the 
following form: 



The first member of equation (24) then becomes 



(B6) 

According to Integral 195 of reference 30, however. 
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therefore 



Equation (B7) Is equation (24) for u * Uq. 


vz.t> 
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Figure 1. - Field of integration for evaluating velocity potential 

(equation (10)) of a source* 



Figure 2. - Relation between time delays T a and Tp and 
position of wave front. 
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Figure 4. - Fields of integration for evaluating upwash between 
wing boundary and foremost Mach wave. 










